Method of determining the yaw rate of a target vehicle

ABSTRACT

This disclosure describes a radar system configured to estimate a yaw-rate and an over-the-ground (OTG) velocity of extended targets in real-time based on raw radar detections. This disclosure further describes techniques for determining instantaneous values of lateral velocity, longitudinal velocity, and yaw rate of points of a rigid body in a radar field-of-view (FOV) of the radar system.

TECHNICAL FIELD OF INVENTION

This disclosure generally relates to a radar system suitable for an automated vehicle, and more particularly relates to a system that calculates a yaw-rate of the target. This invention relates to a vehicle on-board method of estimation of planar motion parameters of objects (targets such as another vehicle) detected by a (e.g. Doppler) radar. It has particular application to detected the planar motion parameters of a vehicle by a host vehicle equipped with a radar system.

BACKGROUND OF INVENTION

Aspects of the invention relate to an improved method of determining the instantaneous values of lateral velocity, longitudinal velocity and yaw rate of any point of a rigid body in the radar field-of-view (FOV). Aspects are applicable (but not limited) to the estimation of vehicles for automotive perception systems and can be used in Active Safety, Driver Assistance and Autonomous Driving applications.

When an automated or autonomous host-vehicle is preceded by a target-vehicle traveling forward of the host-vehicle, it is advantageous for the system that controls the operation (e.g. steering, brakes, engine) of the host-vehicle to have knowledge of the yaw-rate of the target vehicle. Knowledge of the yaw-rate of the target vehicle can be useful to, for example, temporarily allow a closer following distance because the target-vehicle is turning out of the travel path of the host-vehicle.

In prior art, no algorithm was reported for an instantaneous estimation of full planar motion of rigid body target objects based on raw detections of a single radar.

SUMMARY OF THE INVENTION

In one aspect is provided a method of determining the yaw rate ({circumflex over (ω)}_(t)) of a target vehicle in a horizontal plane by a host vehicle equipped with a radar system, said radar system including a radar sensor unit adapted to receive signals emitted from said host vehicle by said target, comprising: emitting a radar signal at a single time-point instance and determining from a plurality (m) of point radar detections measurements captured from said target vehicle by said radar sensor unit in said single radar measurement instance, the values for each point detection of range, azimuth and range rate; [r_(i), θ_(i), {dot over (r)}_(i)]; determining the values of the longitudinal and lateral components of the range rate equation of the target (c_(t), s_(t)) from the results ({dot over (r)}_(i), θ_(i),) of step a); where the range rate equation is

${{\overset{.}{r}}_{i,{cmp}} = {\begin{bmatrix} {\cos\theta_{i}} & {\sin\theta_{i}} \end{bmatrix}\begin{bmatrix} c_{t} \\ s_{t} \end{bmatrix}}},$ where {dot over (r)}_(i,cmp)={dot over (r)}_(i)+u_(s) cos θ_(i)+v_(s) sin θ_(i) where u_(s) is the sensor unit or host vehicle longitudinal velocity and v_(s) is the sensor unit or host vehicle lateral velocity; determining the orientation angle of the target (γ_(t,scs)); d) determining the target center (x_(t) and y_(t)) from the results (r_(i), θ_(i)) of step a); e) determining a line l_(PH) perpendicular to the orientation of the target and passing through the center of the target (x_(t,c,scs), y_(t,c,scs)) from the results of step (c) and (d); f) determining a line l_(CA) passing through the center of rotation of said target and the position of radar sensor unit of said vehicle; from the steps of b); g) determining the intersection point of the lines l_(CA) and l_(PH) from steps e) and f) being the position of the center of rotation [{circumflex over (x)}_(t,COR), ŷ_(t,COR)] of the target; and h) estimating and yaw rate {circumflex over (ω)}_(t) from the position of the center of rotation found in step g) and components of range rate equation of the target (c_(t) or s_(t)) of step b).

Step c) may comprise determining the orientation angle of the target (γ_(t,scs)); from the values of range and azimuth (r_(i), θ_(i),) of said point detections.

Step c) may comprise determining the orientation angle of the target (y_(t,scs)); from L fit, Hough transform or rotating caliper methodology from the results of step a).

Step b) may include estimating components of range rate equation of the target (ct,st) said point detection measurements of azimuth and range rate [θ_(i), {dot over (r)}_(i)].

Said estimating may comprise applying a Least Square method.

Said line l_(PH) of step e) may be determined and defined from the following equation:

$y = {{{{- \frac{1}{\tan\gamma_{t,{scs}}}}x} + {\frac{1}{\tan\gamma_{t,{scs}}}x_{t,c,{scs}}} + y_{t,c,{scs}}} = {{a_{1}x} + {a_{0}.}}}$

Said line l_(CA) of step f) may be determined and defined from the following equation:

$y = {{{- \frac{{\overset{\sim}{s}}_{t,{scs}}}{{\overset{\sim}{c}}_{t,{scs}}}}x} = {{b_{1}x} + {b_{0}.}}}$

In step g) the position of the center of rotation [{circumflex over (x)}_(t,COR), ŷ_(t,COR)] of the target may be determined from the following equations:

${{\hat{x}}_{t,{COR},{scs}} = {\frac{b_{0} - a_{0}}{a_{1} - b_{1}}{and}}},{{\hat{y}}_{t,{COR},{scs}} = {\frac{{a_{1}b_{0}} - {a_{0}b_{1}}}{a_{1} - b_{1}}.}}$

In step h) the yaw rate may be determined from the following equations:

${\hat{\omega}}_{t} = {{\frac{{\overset{\sim}{s}}_{t,{scs}}}{- {\hat{x}}_{t,{COR},{scs}}}{or}{\hat{\omega}}_{t}} = {\frac{{\overset{\sim}{c}}_{t,{scs}}}{{\hat{y}}_{t,{COR},{scs}}}.}}$

The method may include additionally determining estimates of longitudinal velocity û_(t,i), lateral velocity {circumflex over (v)}_(t,i) of certain target point from the value of yaw rate and the co-ordinates of the center of rotation of the target (y_(t,COR,wcs) x_(t,COR,wcs)) using the following equation:

$\begin{bmatrix} {\hat{u}}_{t,i,{scs}} \\ {\hat{v}}_{t,i,{scs}} \end{bmatrix} = {\begin{bmatrix} {\left( {{\hat{y}}_{t,{COR},{scs}} - y_{t,i,{scs}}} \right){\hat{\omega}}_{t}} \\ {\left( {x_{t,i,{scs}} - {\hat{x}}_{t,{COR},{scs}}} \right){\hat{\omega}}_{t}} \end{bmatrix}.}$

BRIEF DESCRIPTION OF DRAWINGS

The present invention is now described by way of example with reference to the accompanying drawings in which:

FIG. 1 shows a target co-ordinate system;

FIG. 2 shows a vehicle coordinate system;

FIG. 3 shows a sensor coordinate system;

FIG. 4 illustrates how to calculate velocity vectors at the locations of three raw detections;

FIG. 5 which shows the result of the cloud algorithm for yawing targets;

FIG. 6 shows an example of the determination of target orientation based on spatial distribution of raw detections;

FIG. 7 illustrates a geometrical interpretation of one method;

FIG. 8 shows an example of test results of a maneuver;

FIGS. 9 a to d illustrate the movement relating to FIG. 8 ;

FIG. 10 shows an example of test results of a maneuver;

FIGS. 11 a to d illustrate the movement relating to FIG. 10 ;

FIG. 12 shows an example of test results of a maneuver; and

FIGS. 13 a to d illustrate the movement relating to FIG. 12 .

DETAILED DESCRIPTION

Reference will now be made in detail to embodiments, examples of which are illustrated in the accompanying drawings. In the following detailed description, numerous specific details are set forth in order to provide a thorough understanding of the various described embodiments. However, it will be apparent to one of ordinary skill in the art that the various described embodiments may be practiced without these specific details. In other instances, well-known methods, procedures, components, circuits, and networks have not been described in detail so as not to unnecessarily obscure aspects of the embodiments.

‘One or more’ includes a function being performed by one element, a function being performed by more than one element, e.g., in a distributed fashion, several functions being performed by one element, several functions being performed by several elements, or any combination of the above.

It will also be understood that, although the terms first, second, etc. are, in some instances, used herein to describe various elements, these elements should not be limited by these terms. These terms are only used to distinguish one element from another. For example, a first contact could be termed a second contact, and, similarly, a second contact could be termed a first contact, without departing from the scope of the various described embodiments. The first contact and the second contact are both contacts, but they are not the same contact.

The terminology used in the description of the various described embodiments herein is for describing embodiments only and is not intended to be limiting. As used in the description of the various described embodiments and the appended claims, the singular forms “a”, “an” and “the” are intended to include the plural forms as well, unless the context clearly indicates otherwise. It will also be understood that the term “and/or” as used herein refers to and encompasses all possible combinations of one or more of the associated listed items. It will be further understood that the terms “includes,” “including,” “comprises,” and/or “comprising,” when used in this specification, specify the presence of stated features, integers, steps, operations, elements, and/or components, but do not preclude the presence or addition of one or more other features, integers, steps, operations, elements, components, and/or groups thereof.

As used herein, the term “if” is, optionally, construed to mean “when” or “upon” or “in response to determining” or “in response to detecting,” depending on the context. Similarly, the phrase “if it is determined” or “if [a stated condition or event] is detected” is, optionally, construed to mean “upon determining” or “in response to determining” or “upon detecting [the stated condition or event]” or “in response to detecting [the stated condition or event],” depending on the context.

Accurate estimation of the yaw-rate and over-the-ground (OTG) velocity is important for many driving-assistance systems. Described herein is a radar system configured to estimate the yaw-rate and OTG velocity of extended targets (largely, for vehicle tracking) in real-time based on raw radar detections (i.e., range, range-rate, and azimuth). As used herein, the term ‘extended-targets’ is used to refer to targets that present multiple, spaced-apart scattering-points so the term ‘extended-target’ is understood to mean that the target has some physical size. The various scattering-points are not necessarily individually tracked from one radar scan to the next, so the number of scatter-points can be a different quantity and/or each scattering point have a different location on the extended-target in successive radar scans.

The invention determines instantaneous values of lateral velocity, longitudinal velocity and yaw rate of any point of a rigid body (such as another vehicle) in the radar field-of-view (FOV). Generally, a host vehicle is equipped with a radar system where reflected radar signals (detection) from another vehicle in the field of view are processed to provide data in order to ascertain these parameters. In order to do this various conditions and requirements are needed. The target (rigid body/vehicle) needs to be a distributed target, i.e. provide a plurality of detections from the same target.

Also assumed is an approximation of the distributed target by a rigid body model which is e.g. appropriate for vehicles (passenger cars, trucks, motorbikes, trains, trams, etc.), though not generally applicable to vulnerable road users.

Radar detections received by the host vehicle form the target provide raw data provide data with respect to the position of the radar transmit/receive element/unit e.g. the Cartesian position of the detections or the Polar coordinates (azimuth angle, range). By using e.g. Doppler techniques, the range rate can also be determined.

In the subsequent concept description, the following conventions and definitions are used:

World Coordinate System

As is convention an inertial coordinate system with the origin fixed to a point in space is used—it is assumed the co-ordinate system does not move and does not rotate. Conventionally the coordinate system is right-handed; the Y-axis orthogonal to the X-axis, pointing to the right; the Z-axis pointing into the page and positive rotation is to the right of the X-axis; see FIG. 1 which shows such a co-ordinate system with origin 1 and a non-ego vehicle (target) 2.

Vehicle Coordinate System

The origin may be located at the center of the front bumper 3 of the host vehicle 4 as shown by FIG. 2 . The X-axis is parallel to the longitudinal axis of the vehicle. The coordinate system is right-handed with the Y-axis orthogonal to the X-axis, pointing to the right, the Z-axis pointing into the page and positive rotation to the right of the X-axis

Sensor Coordinate System

Origin located at the center of the sensor unit/radome. The X-axis is perpendicular to the sensor radome, pointing away from the radome. The coordinate system is right-handed: Y-axis orthogonal to the X-axis, pointing to the right; Z-axis pointing into the page; Positive rotation to the right of the X-axis. FIG. 3 shows a sensor origin 5.

In aspects of the invention and with prior art techniques, the velocity and the yaw rate of the host vehicle is assumed known. The host over the ground (OTG) velocity vector is defined as: V_(h)=[u_(h)v_(h)]^(T), where u_(h)—host longitudinal velocity and v_(h)—host lateral velocity.

Sensor mounting position and boresight angle in the vehicle coordinate system are also assumed known; the following notations are used: x_(s,VCS)—sensor mounting position, longitudinal coordinate; y_(s,VCS)—sensor mounting position, lateral coordinate; and γ_(s,VCS)—sensor boresight angle.

The sensor(s) Over the Ground (OTG) velocities are assumed known (determined from host vehicle motion and sensor mounting positions). Sensor velocity vector is defined as V_(s)=[u_(s) v_(s)]^(T) with u_(s)—sensor longitudinal velocity and v_(s)—sensor lateral velocity. At each radar measurement instance, the radar unit/sensor captures m raw detections from the target. Each raw detection is described by the following parameters expressed in the sensor coordinate system: r_(i)—range (or radial distance), θ_(i)—azimuth angle, and {dot over (r)}_(i)—raw range rate (or radial velocity) i=1, . . . , m.

Target planar motion is described by the Target over-the-ground velocity vector at the location of each raw detection: V_(t,i)=[u_(t,i) v_(t,i)]T where: u_(t,i)—longitudinal velocity at the location of i-th raw detection, and v_(t,i)—lateral velocity at the location of i-th raw detection.

Target planar motion can be described as well by: V_(t,COR)=[ω_(t) x_(t,COR) y_(t,COR)]^(T), where ω_(t)—target yaw rate, x_(t,COR)—longitudinal coordinate of the center of target's rotation, y_(t,COR)—lateral coordinate of the center of target's rotation, x_(t,c)—longitudinal coordinate of the center of target's bounding box, y_(t,c)—lateral coordinate of the center of target's bounding box, and γ_(t)—orientation of the target object

FIG. 4 illustrates how to calculate velocity vectors at the locations of three raw detections (depicted by reference numeral 6) captured from the same rigid body target and the yaw rate of that target. Center of the target's rotation depicted with reference numeral 7. The range rate equation for a single raw detection is given as follows: {dot over (r)}_(i)+u_(s) cos θ_(i)+v_(s) sin θ_(i)=u_(t,i) cos θ_(i)+v_(t,i) sin θ_(i).

To simplify the notation, the notion of a compensated range rate is introduced and defined as: {dot over (r)}_(i,cmp)={dot over (r)}_(i)+u_(s) cos θ_(i)+v_(s) sin θ_(i), with {dot over (r)}_(i,cmp)—range rate compensated of i-th raw detection.

Then the equation is reduced to: {dot over (r)}_(i,cmp)=U_(t,j) cos θ_(i)+v_(t,i) sin θ_(i).

Range rate equation in vector form

${\overset{.}{r}}_{i,{cmp}} = {{\begin{bmatrix} {\cos\theta_{i}} & {\sin\theta_{i}} \end{bmatrix}\begin{bmatrix} u_{t,i} \\ v_{t,i} \end{bmatrix}}.}$

Range rate equation in general vector form

${{\overset{.}{r}}_{i,{cmp}} = {\begin{bmatrix} {\cos\theta_{i}} & {\sin\theta_{i}} \end{bmatrix}\begin{bmatrix} c_{t} \\ s_{t} \end{bmatrix}}},$ with: c_(t)—longitudinal range rate equation solution, and s_(t)—lateral range rate equation solution.

Velocity profile is used as range rate equation solution synonym. Note estimated values are denoted with a hat and Least Square solutions are denoted with a tilde.

Problem Formulation

The problem to be solved can be phrased as follows: estimate velocity vectors {circumflex over (V)}_(t,i) and yaw rate {circumflex over (ω)}_(t) using raw detection measurements [r_(i), θ_(i), {dot over (r)}_(i)] captured from a rigid body target in a single radar measurement instance. Since the locations of the three raw detections are known (by means of direct measurements and sensor mounting position), the equivalent problem formulation is: estimate the position of the center of rotation [{circumflex over (x)}_(t,COR), ŷ_(t,COR)] and yaw rate {circumflex over (ω)}_(t) using raw detection measurements [r_(i), θ_(i), {dot over (r)}_(i)] captured from a rigid body target in a single radar measurement instance.

Cloud Algorithm

Previously the case of a straight-line moving distributed target was considered. This restriction simplifies the estimation problem as the velocity vectors at the location of each raw detections are identical, i.e.: V_(t,i)=[u_(t,i) v_(t,i)]T=[u_(t) v_(t)]^(T)=V_(t), for i=1, . . . , m.

The Cloud Algorithm (CA) was proposed to estimate over-the-ground lateral v_(t) and longitudinal u_(t) velocity of the “cloud” of detections coming from the same target. This was achieved by Least Square solution to the problem defined as follows:

${\overset{.}{r}}_{i,} = {{\begin{bmatrix} {\cos\theta_{i}} & {\sin\theta_{i}} \end{bmatrix}\begin{bmatrix} u_{t} \\ v_{t} \end{bmatrix}}.}$

The algorithm proved to be a very effective technique for instantaneous estimation of target velocity. Additionally, under the additional assumption of a negligible sideslip, the angle {circumflex over (γ)}_(t)=tan⁻¹({tilde over (v)}_(t)/ũ_(t)) can be used as an estimate of the targets heading.

In D. Kellner, M. Barjenbruch, K. Dietmayer, J. Klappstein, and J. Dickmann, “Instantaneous lateral velocity estimation of a vehicle using Doppler radar,” in Proceedings of 16th International Conference on Information Fusion, Istanbul, Turkey, 2013. The same problem and the same theoretical basis for the estimation of lateral velocity of a straight line moving object was considered. Enhancement to the Cloud Algorithm was made by means of executing RANSAC algorithm to identify outliers and executing orthogonal distance regression (ODR) to solve error-in-variables problem for the modified formulation of the original problem. This approach improved robustness of the solution in comparison to the original Cloud Algorithm solution. Computational complexity and the requirement to solve an optimization problem are the major drawbacks of the proposed approach, especially when an application in a production embedded system is to be considered.

Cloud Algorithm Solution for Yawing Targets

Previously the application of the cloud algorithm to the estimation of target's motion without the restriction on straight-line path was investigated. Such situation is shown in FIG. 4 . The over-the-ground velocity vectors at the location of each detection are determined as follows:

$\begin{bmatrix} u_{t,i,{wcs}} \\ v_{t,i,{wcs}} \end{bmatrix} = {\begin{bmatrix} {\left( {y_{t,{COR},{wcs}} - y_{t,i,{wcs}}} \right)\omega_{t}} \\ {\left( {x_{t,i,{wcs}} - x_{t,{COR},{wcs}}} \right)\omega_{t}} \end{bmatrix}.}$

The range rate equation for each raw detection was derived to be: {dot over (r)}_(i,cmp)=(y_(t,COR,wcs)−y_(s,wcs))ω_(t) cos θ_(i)+(x_(s,wcs)−x_(t,COR,wcs))ω_(t) sin θ_(i).

It was then shown that the Least Square solution to this problem results in: ũ_(t)=(y_(t,COR,wcs)−y_(s,wcs))ω_(t), and {tilde over (v)}_(t)=(x_(s,wcs)−x_(t,COR,wcs))ω_(t). One therefore achieves: ũ_(t)(x_(s,wcs)−x_(t,COR,wcs))={tilde over (v)}_(t)(y_(t,COR,wcs)−y_(s,wcs))

This is a result of a major importance. One can conclude (conclusion 1) that by considering range rates and azimuths from a cloud of detections coming from the same rigid body target and captured in a single look of a single radar, it is possible to estimate the position of the center of rotation and the yaw rate. Further (conclusion 2) it is possible to determine the line at which the center of rotation lies. This line passes through the sensor. This is shown in FIG. 5 which shows the result of the Cloud Algorithm for yawing targets. Calculations based on the three raw detections denoted 6 results in the line (of the center of detection) 8 passing through the sensor and the center of rotation 7. The exact location of the center of rotation 7 not specified by Cloud Algorithms.

In the results highlighted above are followed by the proposition of estimating the full motion of the target using detections captured in multiple looks (instances) of a single radar. Such an estimation is possible under the assumption of a constant yaw rate during the time of these multiple looks. It is also necessary for the host vehicle to be moving. Geometrical interpretation of this method is that one needs to find a (stationary) point of intersection of two lines passing through the (moving) sensor in multiple time instances. Because angles of the two lines are likely to be very similar (the faster the host moves the greater the angle), the approach is prone to inaccuracies.

In D. Kellner, M. Barjenbruch, J. Klappstein, jurgen Dickmann, and K. Dietmayer, “Instantaneous full-motion estimation of arbitrary objects using dual doppler radar,” in Proceedings of Intelligent Vehicles Symposium (IV), Dearborn, Mich., USA, 2014, the same problem was considered and the solution to the Conclusion 1 was to take measurements from two sensors. The Conclusion 2 was not mentioned by the authors, but one supposes e that the lines passing through both sensors would intersect in the center of target's rotation. The authors then applied RANSAC and ODR techniques previously presented in order to improve the robustness of the solution.

The drawbacks of the solution proposed above include that the target needs to be in the field-of-view of multiple sensors (in automotive radar configurations, common parts of field-of-view of multiple radars are usually small). For good accuracy and robustness of the solution a significant distance between the two sensors is desired. The scope of such installation is significantly limited by the dimensions of the host vehicle. One either needs to assume constant yaw rate between the looks of two radars or to synchronize the measurement instance of both sensors. Additional data processing applied to make the solution more robust are characterized by significant computational complexity.

Yaw Cloud Algorithm

The concepts reviewed so far have not assumed any motion model of the target. An Extended Kalman Filter has been formulated based on a constant turn motion model and measurement model derived from a cloud algorithm. The algorithm is called YCA (Yaw Cloud Algorithm). In the proposed formulation measurements from a single radar are considered, but the possibility of extension to multiple sensors is also highlighted.

One object in the current invention in comparison to the above, is not to use time-filtering, and not assume any underlying motion model for the target. The approach does not suffer from difficulties associated with appropriate track initialization.

Object Orientation Estimation

Further prior art relevant to the current concept focuses on the estimation of target orientation based on high-resolution radar, laser scanner or LiDAR measurements. These algorithms ignore the range rate measurement and focus exclusively on the locations of raw detections. FIG. 6 shows an example of the determination of target 9 orientation based on spatial distribution of raw detections. Example with three detections denoted with dots 6.

Several algorithms were proposed for such applications. They can be grouped in three main categories: Rotating Calipers based approaches; L-fit; and the Hough transform. Additional modifications to the rotating calipers method were proposed. The accuracies reported for these algorithms based on the assessment of experimental results were about 10 degree RMS. The accuracies were shown to be dependent on the relative target orientation in SCS (sometimes called aspect angle or exposition angle). A further simplified approach suitable for real-time executions has already been implemented for Short Range Radars. The enhancement proposed in that implementation is to treat the orientation of the tracked object associated with a considered cloud as an initial guess for instantaneous orientation estimation. Next the adjustment of the orientation is performed using the distribution of detections captured only in the current look. This approach can be treated as an optional ‘hot starting’ of the instantaneous orientation estimation. The results are promising in particular when the target is close to the host and a sufficient number of raw detections is available. It is noted here that the determination of target orientation in prior art is not related to the estimation of target motion. Two branches of building blocks were identified: development of algorithms considering only azimuth and range rate measurements and development of algorithms considering only azimuth and range measurements

DESCRIPTION OF THE PREFERRED EMBODIMENTS

The invention includes provides instantaneous estimation of target planar motion from detections of a single radar look of e.g. a high resolution Doppler radar with the possibility of detecting that the target has a yaw rate above a certain threshold. It is to be noted that the “raw data” from this single radar look provides the parameters of r_(i)—range (or radial distance), θ_(i)—azimuth angle, {dot over (r)}_(i)—raw range rate (or radial velocity) for each ith point of m point detections on a rigid body. These are the parameters which are used to determine the yaw rate as well as longitudinal and lateral velocities, where i=1, . . . , m. It is to be noted that the term instantaneous or single look radar data would include reflection data from a “chirp” in Doppler techniques which may scan over e.g. up to 2 ms. By this known methodology range rate may be determined.

Methodology according to aspects are based on the combination of the cloud algorithm and spatial (or geometry) based estimation of object heading. In this way all three radar measurements (i.e. range, azimuth and range rate) are exploited to overcome the previously known limitations.

For the sake of convenience, the earlier presented cloud algorithm for yawing targets is now expressed in sensor coordinate system. In an initial step, the method comprises emitting a radar signal at a single time-point instance and determining from a plurality (m) of point radar detections measurements therefrom captured from said radar sensor unit in a said single radar measurement instance the values for each point detection of range, azimuth and range rate; [r_(i), θ_(i), {dot over (r)}_(i)]. In the next step, estimates of range rate equation parameters (of the target) obtained by the Least Square Solution to:

${\overset{.}{r}}_{i,{cmp}} = {{\begin{bmatrix} {\cos\theta_{i}} & {\sin\theta_{i}} \end{bmatrix}\begin{bmatrix} c_{t} \\ s_{t} \end{bmatrix}}.}$

Starting with the definition of over-the-ground velocity vectors at the location of each detection with longitudinal u_(t,i,scs) and lateral velocities v_(t,i,scs) aligned along

${{sensor}{axes}{:\begin{bmatrix} u_{t,i,{scs}} \\ v_{t,i,{scs}} \end{bmatrix}}} = {\begin{bmatrix} {\left( {y_{t,{COR},{scs}} - y_{t,i,{scs}}} \right)\omega_{t}} \\ {\left( {x_{t,i,{scs}} - x_{t,{COR},{scs}}} \right)\omega_{t}} \end{bmatrix}.}$

The range rate equation becomes: {dot over (r)}_(i,cmp)=(y_(t,COR,scs)−y_(t,i,scs))ω_(t) cos θ_(i)+(x_(t,i,scs)−x_(t,COR,scs))ω_(t) sin θ_(i).

Since: y_(t,i,scs) cos θ_(i)=r_(t,i) sin θ_(i) cos θ_(i)=x_(t,i,scs) sin θ_(i) one can write: {dot over (r)}_(i,cmp)=(y_(t,COR,scs))ω_(t) cos θ_(i)+(−x_(t,COR,scs))ω_(t) sin θ_(i). The Least Square solution to this problem is: {tilde over (c)}_(t,scs)=y_(t,COR,scs)ω_(t), {tilde over (s)}_(t,scs)=x_(t,COR,scs)ω_(t), and {tilde over (c)}_(t,scs)(y_(t,COR,scs))={tilde over (s)}_(t,scs)(−x_(t,COR,scs)) It is to be noted that the range measurement is not used in velocity estimation.

In the next step the orientation angle of the target (γ_(t,scs)) must be computed. This may be estimated from raw detection data such as (r_(i), θ_(i),) and may be estimated by L fit, Hough transform or rotating caliper methodology

In the next step the line l_(PH) perpendicular to the orientation of the target and passing through the center of the target (x_(t,c,scs), y_(t,c,scs)) is determined by:

${l_{PH}:y} = {{{{- \frac{1}{\tan\gamma_{t,{scs}}}}x} + {\frac{1}{\tan\gamma_{t,{scs}}}x_{t,c,{scs}}} + y_{t,c,{scs}}} = {{a_{1}x} + {a_{0}.}}}$

In the next step, the line l_(CA) (equivalent to reference numeral 7 of FIG. 5 ) passing through the center of rotation and sensor mounting position is determined by applying the cloud algorithm:

${{l_{CA}:y} = {{{- \frac{{\overset{\sim}{s}}_{t,{scs}}}{{\overset{\sim}{c}}_{t,{scs}}}}x} = {{b_{1}x} + b_{0}}}},$ where b₀ can be zero—so the term may be removed in another form of the equation.

In a fourth step the intersection point of the lines from the previous two steps are determined with:

${{\hat{x}}_{t,{COR},{scs}} = \frac{b_{0} - a_{0}}{a_{1} - b_{1}}},{{{and}{\hat{y}}_{t,{COR},{scs}}} = {\frac{{a_{1}b_{0}} - {a_{0}b_{1}}}{a_{1} - b_{1}}.}}$

It is to be noted that if both lines are parallel (i.e. a₁=b₁), then there is no solution and the yaw rate cannot be estimated. Additionally: if a₁≠b₀: the lines are parallel, there is no intersection thus no yaw rate, the target is moving along a straight line and where a₀=b₀: the lines are identical and nothing can be said about the motion of the target.

In step 5 if both the orientation estimation and cloud algorithm return plausible results, the yaw rate of the target calculate is calculated as:

${\hat{\omega}}_{t} = {{\frac{{\overset{\sim}{s}}_{t,{scs}}}{- {\hat{x}}_{t,{COR},{scs}}}{or}{\hat{\omega}}_{t}} = {\frac{{\overset{\sim}{c}}_{t,{scs}}}{{\hat{y}}_{t,{COR},{scs}}}.}}$

It is to be noted that the estimated yaw rate depends on both the Least Square result for cloud algorithm as well as the position of COR. For example, for the same {tilde over (s)}_(t,scs) yaw rate will be higher when the center of rotation is calculated to be near the host vehicle and will be lower when the {circumflex over (x)}_(t,COR,scs) is further away from the host. FIG. 7 illustrates a geometrical interpretation of the method. Center of rotation estimated to be in the location of the dot 10.

Results

Results of estimating the yaw rate by the proposed technique was evaluated using experimental data. The experiments were collected for stationary host vehicle and a single target vehicle. Both vehicles were equipped with a differential GPS which accurately measures positions, velocities and yaw rates. In post-processing stage, all detections captured from the vicinity of the reference object were selected for calculating cloud algorithm. The orientation of the target was taken directly from reference data.

FIG. 8 shows an example of test results of a maneuver which comprised the target 9 driving on a circular track 12. FIG. 9 a to d illustrate the movement. The point of intersection of the two lines 13 and 14 (equivalent to lines lCA and lPH respectively) remains roughly at the same position, which is consistent with the actual maneuver performed. The only exception is when the two lines are almost parallel. This confirms analytical results.

FIG. 10 shows an example of test results of a maneuver where the target 9 is driving straight and FIG. 11 a to d illustrates temporal movement. The two lines 13 and 14 remain parallel. This is consistent with the expected behavior, because the target is not yawing.

FIG. 12 shows examples of test results of a maneuver where the target is maneuvering at a junction and FIGS. 13 a to d illustrate temporal movement. The intersection point changes its position as the yaw rate of the target smoothly changes during the maneuver. This is consistent with expected behavior. Yet again, promising result for junction assist applications.

In the scenario discussed in this document, detections coming from a single sensor are considered. There are different ways detections from multiple sensors can be used when the target is in the common area of their FOVs: each sensor may estimate its own lines resulting from cloud algorithm and orientation estimation

In ideal case all lines should intersect at the same point, i.e. the center of rotation. Each sensor can calculate its own center of rotation. These points can then be used to determine resulting estimate for the center of rotation. Such processing should consider plausibility measures for sensor estimate. Detections from all sensors can be collected for improved orientation estimation, but each sensor can calculate a cloud algorithm result based on its own detections.

Using Filtered Heading

Estimation of object orientation usually require more detections than the execution of cloud algorithm. For the purpose of yaw rate estimation/detections, one may consider using target orientation from an associated tracker object. In such a case, yaw rate can be estimated based on only two raw detections. The yaw rate estimation is likely to be more robust, but it will not be an instantaneous result.

Include the Proposed Approach in Object Tracking Filter

If instantaneous estimation is not required, one can consider formulating an object tracking filter. The measurement model of such a filter would then include estimations of yaw rate and center of rotation as proposed in this work.

Propagation of Uncertainty

If the uncertainty of cloud algorithm, as well as uncertainty of orientation is known then uncertainty propagation rules can be used for uncertainty estimation of center of rotation and yaw rate. It can be the most efficient way to determinate if motion estimation is reliable. Estimated uncertainty can be also used in Kalman filter as a level of measurement noise. In its pure form, the approach described here is not based on tracking or filtering and therefore does not require any motion model and no initialization stage is required.

The approach has a lower number of restrictions in comparison to currently known concepts based on Cloud Algorithm. It does not need measurements from several time instances and works for stationary host vehicle. It does not need the target to be in the field of view of two sensors mounted significant distance away from each other. The approach does not require filtering/tracking and does not rely on any target motion model. The algorithm is suitable for application in production embedded systems because of its low computational complexity. The method can be immediately used in state-of-art short range radars for estimation of objects at low range (up to about 15 meters); e.g. it is suitable for junction assist applications (crossroads, roundabouts). The accuracy of the results can only be improved if more detections from moving objects were available; this is likely to happen as sensing technology progresses and in sensors dedicated to Autonomous Driving applications. 

We claim:
 1. A method of determining a yaw rate ({circumflex over (ω)}_(t)) of a target vehicle in a horizontal plane by a host vehicle equipped with a radar system, said radar system including a radar sensor unit adapted to receive signals emitted from said host vehicle and reflected by said target vehicle, estimated values being denoted by a hat ({circumflex over ( )}) said method comprising: a) emitting a radar signal at a single time-point instance and determining, from a plurality (m) of point radar detections measurements captured from said target vehicle by said radar sensor unit in said single time-point instance, values for each point radar detection including range, azimuth, and range rate: [r_(i), θ_(i){dot over (r)}_(i)], to provide results of step (a); b) determining, by using a Least Square method and in a sensor coordinate system, Least Square solution values of longitudinal and lateral components of a range rate equation of the target vehicle ({tilde over (c)}_(t,scs),{tilde over (s)}_(t,scs)) from the results ({dot over (r)}_(i),θ_(i),) of step a), where the range rate equation is ${\overset{.}{r}}_{i,{cmp}} = {\begin{bmatrix} {\cos\theta_{i}} & {\sin\theta_{i}} \end{bmatrix}\begin{bmatrix} c_{t} \\ s_{t} \end{bmatrix}}$ and {dot over (r)}_(i,cmp)={dot over (r)}₁+u_(s) cos θ_(i)+v_(s) sin θ_(i), u_(s) is the radar sensor unit's velocity or the host vehicle's longitudinal velocity and vs is the radar sensor unit's velocity or the host vehicle's lateral velocity, to provide results of step (b), an x-axis of the sensor coordinate system being perpendicular to the radar sensor unit and pointing away and ahead of the radar sensor unit, a y-axis of the sensor coordinate system being perpendicular to the x-axis and pointing right of the radar sensor unit, Least Square solution values being denoted by a tilde ({tilde over ( )}) and solution values in the sensor coordinate system being denoted with a scs subscript (_(scs)); c) determining an orientation angle (γ_(t,scs)) of the target vehicle relative to the x-axis of the sensor coordinate system and another line representing an orientation of the target vehicle, to provide results of step (c); d) determining a center of the target vehicle (x_(t,c,scs) and y_(t,c,scs)) from the range and azimuth (r_(i), θ_(i)) of step a), to provide results of step (d); e) determining a line l_(PH) perpendicular to the other line representing the orientation of the target vehicle and passing through the center of the target vehicle (x_(t,c,scs) and y_(t,c,scs)), the determining from the results of step (c) and (d), respectively, to provide results of step (e); f) determining a line l_(CA) that passes through a center of rotation of said target vehicle using a position of the radar sensor unit of said host vehicle and the longitudinal and lateral components of the range rate equation of the target vehicle ({tilde over (c)}_(t,scs),{tilde over (s)}_(t,scs)), from the result of step b), to provide results of step (f); g) responsive to a determination that the line l_(CA) and l_(PH) are not parallel, determining an intersection point of the lines l_(CA) and l_(PH) from steps e) and f), the intersection point being a position of the center of rotation [{circumflex over (x)}_(t,COR,scs), ŷ_(t,COR,scs)] of the target vehicle, to provide results of step (g); h) estimating the yaw rate ({circumflex over (ω)}_(t)) based on the position of the center of rotation found in step g) and based on the components of the range rate equation of the target vehicle ({tilde over (c)}_(t,scs),{tilde over (s)}_(t,scs)) of the results of step b); and i) operating, by a driving-assistance system of the host vehicle, the host vehicle based on the yaw rate of the target vehicle.
 2. A method as claimed in claim 1, wherein said step c) comprises determining the orientation angle (γ_(t,scs)) of the target vehicle from the values of range and azimuth (r_(i), θ_(i)) of said point detections.
 3. A method as claimed in claim 1, wherein said step c) comprises determining the orientation angle (γ_(t,scs)) of the target vehicle from L fit or Hough transform.
 4. A method as claimed in claim 1, wherein said line l_(pH) of step e) is determined and defined from the following equation: $y = {{{{- \frac{1}{\tan\gamma_{t,{scs}}}}x} + {\frac{1}{\tan\gamma_{t,{scs}}}x_{t,c,{scs}}} + y_{t,c,{scs}}} = {{a_{1}x} + {a_{0}.}}}$ where $a_{1} = {- \frac{1}{\tan\gamma_{t,{scs}}}}$ and $a_{0} = {{\frac{1}{\tan\gamma_{t,{scs}}}x_{t,c,{scs}}} + {y_{t,c,{scs}}.}}$
 5. A method as claimed in claim 4, wherein said line l_(CA) of step f) is determined by applying a cloud algorithm, the cloud algorithm defined as: $y = {{{- \frac{{\overset{\sim}{s}}_{t,{scs}}}{{\overset{\sim}{c}}_{t,{scs}}}}x} = {{b_{1}x} + b_{0}}}$ where {tilde over (s)}_(t,scs)=−{circumflex over (x)}_(t,COR,scs){circumflex over (ω)}_(t); {tilde over (c)}_(t,scs)=−ŷ_(t,COR,scs){circumflex over (ω)}_(t); ${b_{1} = {- \frac{{\overset{\sim}{s}}_{t,{scs}}}{{\overset{\sim}{c}}_{t,{scs}}}}};$ and b₀=0.
 6. A method as claimed in claim 5, wherein in step g) the position of the center of rotation [{circumflex over (x)}_(t,COR,scs), ŷ_(t,COR,scs)] of the target vehicle is determined from the following equations: ${{\hat{x}}_{t,{COR},{scs}} = {\frac{b_{0} - a_{0}}{a_{1} - b_{1}}{and}}},{= {\frac{{a_{1}b_{0}} - {a_{0}b_{1}}}{a_{1} - b_{1}}.}}$
 7. A method as claimed in claim in 6, wherein in step h) the yaw rate is determined from the following equations: ${\hat{\omega}}_{t} = {{\frac{{\overset{\sim}{s}}_{t,{scs}}}{- {\hat{x}}_{t,{COR},{scs}}}{or}{\hat{\omega}}_{t}} = {\frac{{\overset{\sim}{c}}_{t,{scs}}}{{\hat{y}}_{t,{COR},{scs}}}.}}$
 8. A method as claimed in claim 7 including additionally determining estimates of longitudinal velocity û_(t,i), and lateral velocity {circumflex over (v)}_(t,i) of a certain point of the target vehicle from the value of yaw rate and coordinates of the center of rotation of the target vehicle using the following equation: $\begin{bmatrix} {\hat{u}}_{t,i,{scs}} \\ {\hat{v}}_{t,i,{scs}} \end{bmatrix} = {\begin{bmatrix} {\left( {{\hat{y}}_{t,{COR},{scs}} - y_{t,i,{scs}}} \right){\hat{\omega}}_{t}} \\ {\left( {x_{t,i,{scs}} - {\hat{x}}_{t,{COR},{scs}}} \right){\hat{\omega}}_{t}} \end{bmatrix}.}$
 9. The method as claimed in claim 1, wherein the other line representing the orientation of the target vehicle passes through the center of the target vehicle (x_(t,c,scs) and y_(t,c,scs)). 